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We investigate two interacting open quantum systems whose time evolutions are governed by 
Markovian master equations. We show a class of coupled systems whose interaction leaves invariant 
the steady state of one of the systems, i.e., only one of the reduced steady states is sensitive to the 
interactions. A detailed proof with the help of the Trotter product formula is presented. We apply 
this general statement to a few models, one of which is the optomechanical coupling model where 
an optical cavity is coupled to a small mechanical oscillator. 
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I. INTRODUCTION 

Quantum Markovian master equations are the simplest 
but also the most used cases of study for the dynam¬ 
ics of open quantum systems. The concept of open sys¬ 
tems plays an important role since perfect isolation in 
the quantum realm and access to all degrees of freedom 
is not possible. The main source of inspiration was the 
theory of lasers where models describing the interaction 
of a central subsystem weakly coupled to an uncontrol¬ 
lable environment were studied. In the last two decades 
there has been great interest in the study of compos¬ 
ite quantum ^sterns in order to create entanglement be¬ 
tween them (l| with the ultimate aim to realize quantum 
information procedures 0. However, there is always an 
environment causing decoherence which is also the bot¬ 
tleneck to realize highly efficient quantum information 
protocols. Therefore the theory of open quantum sys¬ 
tems was frequently applied to various models Q. 

There is a subclass of models which can be decomposed 
into many interacting parts which in turn are coupled to 
independent environments. The best example is the case 
of two-level atoms interacting with a single mode of the 
radiation field and coupled to independent reservoirs, like 
the model of the many-atom laser or the one-atom 
maser Q. The time evolution of the composite system is 
derived from a microscopic model or is assembled on the 
ground of phenomenological considerations from equa¬ 
tions which describe independently the parts. The latter 
case usually assumes a low coupling strength between the 
two parties which allows a separate description of dissi¬ 
pative effects. In both cases the possible treatments of 
environments consist in Markovian and non-Markovian 
approaches 0. Obtaining the full solution to the dy¬ 
namical equation is a very complicated task even with a 
Markovian approach [3 and sometimes one has to be 
satisfied by identifying the steady state Q or its symme¬ 
try properties In general, the role of the interaction 
in a composed system was studied typically for two cou¬ 
pled harmonic oscillators [13, [Hi and for two interacting 
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identical atoms (l^ . 

In this paper, we focus on the case when two indepen¬ 
dent quantum systems A and B are subject to Markovian 
master equations and their interaction is studied through 
an equation which is obtained by adding the two master 
equations to a Hamiltonian interaction. For physical sys¬ 
tems, this consideration usually implies that the interac¬ 
tion strength between A and B has to be the smallest 
parameter in the system. We concentrate on a special 
class of interaction Hamiltonians that cause dephasing 
on system A and lead to a generator (Liouvillian) of the 
complete system that commutes with the Liouvillian of 
system A’s master equation. The motivating example is 
the optomechanical coupling between a single-mode elec¬ 
tromagnetic field and a small mechanical oscillator in the 
presence of independent decoherence mechanisms [l3 | . It 
will be demonstrated that the considered class of inter¬ 
action Hamiltonians leads to an invariant steady state of 
system A, i.e., tracing out system B in the steady state of 
the composite system one always obtains the same steady 
state of system A which is independent of the interaction 
Hamiltonian. In contrast, the steady state of system B 
may depend on the interaction. In the case of finite di¬ 
mensional systems, theorems of linear algebra can be em¬ 
ployed to prove the statement. However, when dealing 
with infinite dimensional systems a different and more 
sophisticated method is required. Therefore, the proof 
of the invariant steady state is based on the properties 
of quantum dynamical semigroups where we employ the 
Trotter product formula and identify the growth bounds 
of the master equations. While the proof of our statement 
is rather general we will explicitly show it in some simple 
models. The steady state of system B does not have any 
particular property and its solution is model dependent. 
Therefore in our examples we study the steady state of 
system B by analyzing its average excitation. Our work 
presents a general approach towards composite systems 
governed by independent Markovian master equations. 

This paper is organized as follows. In Sec. [Tllwe discuss 
the motivation of our analysis. In Sec. IHII we show the 
general statement about the steady state. In Sec. IlYl we 
apply these findings to three different systems: two inter¬ 
acting spins, a spin interacting with a harmonic oscillator 
and two coupled harmonic oscillators. A discussion about 
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our findings is summarized up in Sec. |Vl In Appendix 1X1 
we briefly discuss the Liouvillians of the Markovian mas¬ 
ter equations involving some technical problems regard¬ 
ing infinite dimensional Hilbert spaces. In Appendix |B] 
we analyse the form of the equilibrium state of a damped 
harmonic oscillator in the number state representation. 


II. MOTIVATION 


There are many physical situations in which the dy¬ 
namics of the system can be decomposed into two par¬ 
ties subject to a Hamiltonian interaction. The individual 
sub-systems might include a relaxation mechanism that 
under the Markovian assumption leads to a master equa¬ 
tion for the density matrix g of the composite system 
that can be written as 

^ = Cg = C^g +C^g + i[g,Hi]. (1) 

The superoperators and are generators of the in¬ 
dividual systems and Hj is the interaction Hamiltonian. 
One system that allows this type of separation and has 
gained a lot of recent interest is the model of optomechan¬ 
ical coupling, where the radiation held and a mechanical 
oscillator are coupled via radiation pressure. In the sin¬ 
gle mode scenario, numerical simulations show that the 
steady state of the radiation field remains independent 
of the coupling strength and of the parameters of the 
mechanical oscillator. This interesting effect can be gen¬ 
eralized to two systems where one of them is composed 
of spins and harmonic oscillators and the other system 
stays arbitrary and the interaction between them holds 
the same symmetries as the optomechanical coupling. In 
the following we motivate the problem by showing the 
optomechanical model. We consider the simplest Hamil¬ 
tonian describing the interaction between a single mode 
of the radiation held and a vibrational mode of a me¬ 
chanical oscillator. This interaction is derived for a cavity 
with a moving mirror mill and the quantization of the 
radiation field is considered with time-varying boundary 
conditions. The single mode assumption for both quan¬ 
tized held and mirror motion was successfully applied to 
describe most of the experiments to date and leads to the 
following interaction Hamiltonian (fi = I) 

Hi = ga^a{b^ +b). ( 2 ) 

Thereby, a (a^) is the annihilation (creation) operator of 
the single mode radiation held with frequency w and b 
( 6 ^) is the annihilation (creation) operator of the mov¬ 
ing mirror’s mode with frequency v. The parameter g 
describes the optomechanical coupling strength. This in¬ 
teraction Hamiltonian commutes with the free Hamilto¬ 
nian of the radiation held uja^a but it does not commute 
with the free Hamiltonian of the mirror iibH. Account¬ 
ing to an irreversible damping mechanism the Liouville 


operators of the individual systems can be cast in the 
following form 

C^g = iuj[g,a^ci\+K{n+l)'D[a\g + KfiD[a^]g, (3) 
g = iv[g,y^b\ + ^{m+l)V^]g + ^friD^^]g, 

which are written in terms of the Lindbladian dissipator 

Vla]^ = aga^ — ag + ga^ a). (4) 

The constant k is the loss rate of the radiation field and 
7 is the damping rate of the mirror’s mode. Both are 
assumed to be coupled to a thermal reservoir with an 
average number of photons n and phonons m. 

An interesting feature can be distinguished fairly easily 
in this model by evaluating the mean occupation photon 
and phonon numbers at the steady state of the dynamics. 
This calculation is explained in detail in Sec. IIVI and here 
we merely show the result, namely 

(d^a)st = h, (5) 

An^g'^ 4h(h-I-1)(2 k-I- 7 ) 5 ^ 

7 ^ -I- Av'^ 7 [(2k -|- 7)2 + 4i/2] 

The most striking effect is that the mean photonic occu¬ 
pation number remains unchanged, i.e. it is insensitive 
to the coupling to the mechanical oscillator. In contrast 
the mean phonon number grows with increasing photon 
number n and coupling strength g. In the special case 
when the photonic bath is at zero temperature (h = 0 ), 
the occupation numbers of the cavity and the moving 
mirror remain unchanged. It is not a complicated task 
to show that the composite steady state is insensitive to 
the interaction in this situation. This is actually the ap¬ 
proximate situation of current experimental scenarios, as 
for optical frequencies h can be regarded to be zero. In 
addition, the phononic contributions due to the coupling 
in the second line of Eq. ([5]) are negligible in the low cou¬ 
pling regime for low values of h. Nevertheless, with the 
new generations of optomechanical systems reaching the 
strong coupling regime and for microwave frequencies, 
this effect can become important. It is worth mentioning 
that if the coupling is too strong, a more detailed model 
of the loss mechanisms has to be considered such as in 
Ref. [T^ where the eigenstates of the Hamiltonian in Eq. 
© where used to derive the corresponding master equa¬ 
tion. It can be shown that even in this case the effect 
persists: the mean photon number remains unchanged. 
We will comment more on this at the end of Sec. m 
In an attempt to explain the effect one can realize that 

(C^g) e*^^‘ = . (6) 

It turns out that by identifying this special characteris¬ 
tic one is able to generalize this behaviour. There are 
two main purposes of this work. One is to demonstrate 
that for h 7 ^: 0 the steady state of the optical mode re¬ 
mains unchanged in contrast to the changing mechanical 
mode. Second, we will extend the result to a broader 
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class of systems that share the same symmetry relations. 
We will show that these arguments result in an invariant 
steady state for the system with symmetries and this is 
independent of the details of the other system. 

III. STEADY STATE ANALYSIS 

In this section we generalize the statement presented 
in Sec. m and prove it by using the properties of dy¬ 
namical semigroups. First, let us consider two quan¬ 
tum systems A and B, where each of these systems 
undergo a Markovian non-unitary dynamics, described 
by the Gorini-Kossakowski-Sudarshan-Lindblad genera¬ 
tor [iMl- The non-unitary evolution is a trace of the 
interaction of each system with an independent environ¬ 
ment. These environments are considered to be such that 
the evolution of systems A and B has a Markovian irre¬ 
versible behaviour. The dynamical equation of the joint 
system dg/dt = Lg is governed by the Liouville operator 

Cg = C^g + C^g-i[VH^,Q]. (7) 

We assume that system A is composed of Np harmonic 
oscillators and Ns spins and that its dynamics is de¬ 
scribed by the Lindblad master equation 

Np Ns 

= ( 8 ) 

p—l S = 1 

The operator N denotes the reduced density matrix of 
system A alone and the individual Liouvillian for each 
oscillator has the form 

^HO = + KphpT>[a]i], (9) 

which is written in terms of the dissipator defined in Eq. 
Q. dp (ajj) is the annihilation (creation) operator of the 
pth harmonic oscillator. Similarly we have for each spin 

= 7s(l - rhs)'D[d-'"^'^] + (10) 

with the ladder operators ± icF\p)/2 ex¬ 

pressed in terms of the Pauli spin operators for the sth 
spin. All the constants 7 p and Kg have dimension of fre¬ 
quency and are assumed to be non-zero. The dimen¬ 
sionless parameters {hp > 0}p=i_,,_Arp are proportional to 
the temperature of the individual thermal reservoir to 
which each oscillator is coupled and 0 < rhs < 1 for all 
s G {1,..., Ns\. 

It is a challenge in the description of a harmonic os¬ 
cillator that we have to work with unbounded operators, 
which cannot be defined everywhere. This issue is inher¬ 
ited by the Liouvillian in Q which is defined on elements 
of a Banach space. In order to discuss the properties of 
this Liouvillian we require tools from the theory of dy¬ 
namical semigroups. These tools are supported by Ap- 
oendixIAl where we show the spectrum of the Liouvillians 


and we discuss the definition of the adjoint. Examples 
related to the domain of the harmonic oscillator’s Liou¬ 
villian are also presented. 

We assume for system A the following form of the in¬ 
teraction Hamiltonian 

Np Ns 

= ^dldp + ^atN ■ ( 11 ) 

p—l 

The operator Va can be regarded as the excitation op¬ 
erator of system A and therefore it is a dephasing type 
interaction. Our last assumption is that the Hamiltonian 
of system A commutes with the interaction Hamiltonian, 

i.e. 

= ( 12 ) 

Examples of are the Jaynes-Cummings, Ising, and 
Tavis-Cummings Hamiltonians. If these conditions are 
met, we can state that the steady state of the reduced 
system A is independent of the interaction with system 
B. This can be restated as follows: if Pst is the steady 
state of the complete system and solution of the equation 
Cpst = 0, then 

= 0, (13) 

and therefore Nt = Tr^l^st} with NNt = 0- This state¬ 
ment does not hold in general for system B as 

/:^TrA{&t} ^ 0. (14) 

To prove this statement let us start by rewriting the 
Liouvillian £ in terms of superoperators in the following 
form 

L = L^ +N-iC^- - iNA^ (15) 

where 

Cg=lV\gl A^g=^{V\g}, *G{A,i?}, 

C^g=[H^,g]. (16) 

One can note that the superoperator commutes with 

the Liouvillian, 

[C^,£]=0. (17) 

The commutation with the first and last three terms in 

Eq. (ITSl) is evident from Eqs. (fT^ and (ITCl) . and be¬ 
cause superoperators acting separately on systems A and 
B commute. The commutation with N can be shown 
using the commutation relation between the interaction 
Hamiltonian and the Lindblad operators. 

In order to exploit the commutation property between 
£ and C^, it is convenient to work in the eigenbasis of 
the superoperator which is defined in (ITCl) . As it is 
formed by the operator it is convenient to identify 
first the eigenvalue equation of i.e. 

j) = n|n,j), nJeN. (18) 
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Due to © we have 

V-^H^\n,j)=nH^\n,j) (19) 


In search of the steady state structure restricted to sys¬ 
tem A we consider the reduced dynamical equation and 
we note that for the subspace labeled by Z = 0 one has 


which means that does not couple vectors with dif¬ 
ferent n. 

The kets \n,j) represent states where the number of 
photons plus the number of excited spins sum to n. These 
states have a degeneracy d„ = dn{Np,Ns) that is enu¬ 
merated by the index j = and depends on 

the number of spins and oscillators in system A. For in¬ 
stance, in the Jaynes-Cummings model these states are 
|n, 1) = \n) (g) |e) and |n, 2) = |n-|-l) 0 \g). Using the 
states in (flSl) one can find that the eigenvectors of the 
superoperator are given by 


O';,!/ = 


1 ^ 1 +^ o 

2 ’J 


\±d 
2 ’ 


iy = {n,j,k}, (20) 


with I G Z and where we have introduced the collective 
index v. These eigenvectors are rank 1 operators and it 
can be verified that they solve the eigenvalue equation 

= ( 21 ) 

with the eigenvalue Z £ Z. As these eigenvectors are 
degenerate, any superposition of with fixed I is also 
an eigenvector of Furthermore, they are orthogonal 
with respect to the Hilbert-Schmidt inner product 

= TrA^a-j ( 22 ) 

where we have introduced for convenience a ket-bra no¬ 
tation and to write the inner product. Let 

B{'H) be the set of all bounded linear operators on a 
Hilbert space TL. Given the fact that the states \n,j) 
are complete in the Hilbert space it follows that the 
eigenvectors form a complete orthonormal basis for 
the Hilbert space of Hilbert-Schmidt operators 

B 2 {n^) := {x G B{n^) : Tr{iti} < oo} (23) 

which contains the set of all density matrices of system 
A. Therefore, in the previously defined ket-bra notation 
one can write the completeness relation as 

with Vi = (24) 

l V 


VoTveiCg} = [c^ - iPoTrsI^}. (27) 

We have used Eqs. m and (EID, and the relations 
Tr_B{U^p} = TrsIC^p} = 0. This shows that if gst 
is the steady state of C and is the steady state of 
then 


Vopi = VoTrB{gst}. (28) 

This still does not imply the statement in Eq. m as 
we have to show that the equality also holds for the sub¬ 
spaces I ^ 0. 

Eor the remaining part of the proof we use a more 
technical approach based on known results of quantum 
dynamical semigroups [2^. In this sense, let us recall 
that the Liouville operator C in Eq. o is the genera¬ 
tor of the quantum dynamical semigroup . These two 
operators act on density matrices g which are trace-class 
operators in the Hilbert space H = 0 'H®. The Ba¬ 

nach space of all self-adjoint trace-class linear operators 
containing the set of density matrices can be defined as 
0 


m) := {g G Bin) : g^=g, ||^||i < oo}. (29) 
with the norm 


ll^lli =Tr|v^|. (30) 

It was shown by Kossakowski (see Definition I and The¬ 
orem I in [2^) that the Liouville operator £ generates 
a contraction semigroup on Lin). This means that the 
norm of the evolution operator is bounded as ||e'^*||op < 1 
which implies that the norm of any evolved element 
e^^g G Lin) is bounded as 

||e^‘e||l<||e^‘||op||^||l<||^.||l,V^,GL(7^) (31) 

where we introduced the operator norm 

||(!I||op=sup{||(!U;||i: x G Lin), ||x||i = l}. (32) 


Thereby we have introduced the projectors Vi that 
project onto the subspace of fixed value of the eigen¬ 
value I of C^. Using the completeness relation one can 
separate the Liouville operator in different subspaces as 
C = '^I'PiC where we have 

Vi£ = £Vi. (25) 

Vi£ acts only onto the subspace containing eigenelements 
ijpi, of for fixed value of Z. In particular this superop¬ 
erator can be written in terms of the projectors as 

ViC^ = C'^Vi = IVi- (26) 


Our goal now is to find a smaller bound than the bound 
given in Eq. (IHTl) for the elements git) = but pro¬ 
jected to the subspaces with fixed value |Z|. Eor this pur¬ 
pose we consider the projectors 

Qi=Vi+V-u Z^O, (33) 

which act as the identity operator for system B. By 
taking fixed values of |Z| one ensures that the self-adjoint 
property is preserved, i.e. (Qzp)^ = Q;/3 whenever p'l’ = p. 
This does not hold when projecting with Therefore 

Qig G Lin) and the operator norm of each projector 
llQillop = 1- The projector Qi projects all elements x G 
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L{T-L) onto a closed linear subspace oi L{'H) which is also 
a Banach space with norm ||.||i: 


Hamiltonian . This is a stronger result than the 

original requirement for our statement, because 


L\%) :={Qix: x € L{n)} . 


(34) 


TrB{Qie(t -)■ oo)} = 0 TrB{e(t -t oo)} = (41) 


In the next step we take advantage of the fact that 
the generators C and do not mix spaces of different 
value of Z, i.e. Qi£ = CQi- This allows us to define the 
operators restricted to L^{H) 

Cl ■■= QiC = CQu Cf ■.= QiC^ = C^Qi. (35) 

Now we can evaluate the norm of the projected elements 
as 


\\Qimh<\\e^‘%JQim\\i- (36) 

where we also used that Qf = Qi. The operator norm 
11.11 op is induced by the Banach space L^H). In order to 
find a better bound we take advantage of our knowledge 
of the generator C^ and separate the dynamics with the 
use of the Trotter product formula (see Corollary 5.8 in 
Chap. Ill of Ref. (^ 1 




hm (37) 

N^oo V / 


As discussed in Appendix |A] the Liouvillian C^ has 
only a point spectrum . Therefore in the Banach space 
the dynamical semigroup of system A has the op¬ 
erator norm || exp(£)^t)||op < e’*'* (see Definition 2.1 and 
Corollary 3.12 in Chap. IV of Ref. [2l|) where r]i € M_ is 
the largest eigenvalue of Cf in the Banach space L^T-L) 
{I ^ 0). This follows from the fact that there is only one 
zero eigenvalue of C^ with corresponding eigenvector 
that has no projection on the subspaces with I ^ 0, i.e. 
Qip^t = 0 and Vopfi = p^. The rest of the eigenvalues of 
C^ are negative. This allows us to find a bound indepen¬ 
dent of N for the evolution operator that we expanded 
using the Trotter formula in Eq. dSZl)- Using the submul- 
tiplicative property of the operator norm ||.|lop and the 
fact that {Cl — Cf) also generates a strongly continuous 
one-parameter contracting semigroup, we have 

< l|e^^*/^||op < (38) 

which yields 

II lim (39) 

AT—>00 V / '■ 


but we found that Qig{t 00) = 0. 

The above approach can also be used for system A 
alone in the case when we know the spectrum of the Li- 
ouville operator, but due to the internal interactions it 
is hard to solve the dynamics or to find the steady state. 
As an example someone can think about the Jaynes- 
Cummings model where the two-level system and the 
harmonic oscillator have an independent damping mech¬ 
anism. Despite the fact that finding the analytical form 
of the steady state is a challenging task, in an appropri¬ 
ate representation we have to deal only with the block 
diagonal part of the density matrix. 

We end this section with a comment on the choice of 
separate damping mechanisms. This is a valid descrip¬ 
tion for weakly interacting systems, where the environ¬ 
ments are insensitive to the dressed eigenstates of coupled 
systems A and B. For strong coupling, the physical de¬ 
scription is not valid anymore and a more detailed lossy 
model has to be derived. For instance, in Ref. [l^ this 
was investigated for the optomechanical model and the 
master equation was derived using the dressed state of 
the Hamiltonian in Eq. ©• The result is a master equa¬ 
tion that differs from Eq. o and m by an additional dis¬ 
sipative term proportional to D[d^d] and where the me¬ 
chanical dissipators are replaced by D[&] — >■ Vlb — g/i/a)a] 
and 'D\w] —>■ D[U — g/va^a]. With these corrections, the 
condition in Eq. (HB is still fulfilled and therefore the 
steady of system A, i.e. the optical mode, remains invari¬ 
ant as well. In general, our method applies to any two 
dissipative systems A and B with a separable Hamilto¬ 
nian interaction and where: the assumptions (fl^ 

and (HB are satisfied, and one can identify a generator 
C^ acting solely on system A whose zeroth eigenvalue be¬ 
longs only to the eigenspace / = 0 of (see Eq. (HB)- 
There are techniques to infer about the spectrum of a fi¬ 
nite dimensional operator. This implies that our results 
can be applied to finite dimensional systems including 
corrections due to the strong coupling limit. However, 
in the case of infinite dimensional Hilbert spaces with 
non-zero temperature baths, the choice is restricted be¬ 
cause only the spectral decomposition of the optical mas¬ 
ter equation (damped harmonic oscillator) is known. 


IV. EXAMPLES 


This upper bound is independent of N which implies the 
following bound for the norm 

||Qie(Z)l|i < e'"‘||£>(0)||i ^ lim Qig(t) = 0. (40) 

t—¥oo 

Thus the steady state has projection 0 on the subspaces 
I 7^ 0. This fact does not depend on the interaction 


A. A case study of two interacting spins 

The situation can be easily understood for the sim¬ 
plest case of two interacting spin systems {A and B) with 
independent Markovian dynamics. We consider an in¬ 
teraction Hamiltonian which satisfies the symmetries of 
the optomechanical system. Due to the simple structure 
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of the Hilbert spaces the steady state of the whole sys¬ 
tem has an analytical formula. The Hilbert spaces are 
'Ha — 'Hb = with orthonormal basis {|0 )a, |1)a} for 
system A and {|0 )_b, |1)b} for system B. Each system 
alone is governed by the master equation djA/dt = 
with i G {A, B} expressed in terms of the the Liouvillians 

Cp^ = zuj[p\ ai] + Y(1 - s^)V[a^_]p^ + s^V[d\]p\A2) 

where w* is the frequency difference between the energy 
levels | 0 )i and |l)i. 7 ® has the dimension of frequency 
and 0 < s* < 1. The spin operators are given by 

d: = |l).(l|i - |0).(0|., dL = (dV)t = |0).(1|.. (43) 

The steady state of each separate system has the same 
form 

p*t = (l-s*)|0).(0|. + s*|l).(l|.. (44) 

By considering the interaction Hamiltonian 

Hi = (df + d^) , (45) 

one can realize that the interaction Hamiltonian com¬ 
mutes with the free Hamiltonian of system A, i.e. 

lHi,af]=0. (46) 


The steady state pst of the whole system can be analyti¬ 
cally evaluated due to the small dimension of the Hilbert 
space. We do not show the explicit expressions here but 
simply emphasize that it has the form 


^st — 


/ Po.0,0.0 Po.0,0.1 
Po,1.0.0 Po,1.0.1 
0 0 

Vo 0 


0 0 

0 0 

Pi.0.1,0 Pi.0.1,1 
Pi,1,1,0 Pi,1,1,1 


(47) 


where Pij,k,i = {i\A{j\BQst\k) a\1)b and i,j,k,l G {0,1}. 
The reduced steady state of system A is given by 

TrB{est}= ^A^ =pft, (48) 

whereas for system B it reads 

TrA{Pst}= (^§1 ^pff (49) 

The steady state of system B is given by the real pa¬ 
rameter Bi and the complex parameter B 2 both of them 
being functions of w, 7 * and s* {i G {A,B}). In the 
case of strong interaction between systems A and B the 
following relations hold 


limHi = l/2, limH2 = 0, (50) 

^—>•00 Cl—^OC 


resulting in a maximally mixed state for Tr^{pst}- 
In Fig. [T] we found that the population inversion 
(df )st = Trjdf gstj of system B has a monotone decreas¬ 
ing behaviour as a function of the interaction strength H, 


0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 
0.0 

FIG. 1: The population inversion of system A and B in the 
steady state as a function of H/ 7 '^. We set = 1, 

s® = 0.6 and = 0.8. The solid line shows that in system 
A the average excitation difference between the energy levels 
|1)a and |0)a remains independent of interaction strength H, 
while the dashed and the dotted line related to system B 
shows a decrease. The parameters are = 1 (dotted 

curve) and 01 / 7 '^ = 10 (dashed curve). 
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and furthermore its dependence on the frequency differ¬ 
ence between the energy levels uj governs the rate of the 
decrease. The decrease is slowed down by the increase of 
the frequency difference u. System A returns to the same 
steady state independently of the interaction strength fl. 
The steady state of system B depends on the interaction 
strength which can outplay the relaxation mechanism of 
system B. This preliminary result already hints that 
this model can only explain physical situations whereas 

UJ. 


B. A spin system interacting with a harmonic 
oscillator 

In this section we consider an interaction Hamiltonian 
between a finite and an infinite dimensional system. The 
Hilbert spaces are Ha = and Hb = rs(C) (sym¬ 
metric Fock space) with orthonormal basis {|0)a, |1)a} 
for system A and {|n)B}jigN for system B. System A is 
governed by the master equation 

= iwa[p^, -b 7^(1 - s)'D[a-]p^ + 7^sT>[(T+]/3^, 

where we use the spin operators dz = | 1 )a( 1 U — | 0 )a( 0 |a 
and d- = d^ = |0)a(1|a. The parameter wa stands for 
the frequency difference between the levels | 1 )a, | 0 )a, 7 ^ 
has the dimension of frequency and 0 < s < 1. System 
B is governed by the master equation 

^ = C’P” ( 52 ) 

= iujB[p^,b'^b] -b 7'®(1 -b n)'D[b]p^ + 'y^n'D[b^]p^, 
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where b {b^) is the annihilation (creation) operator of 
the harmonic oscillator with frequency ujb- h > 0 is 
proportional to the temperature of the thermal reservoir 
to which the harmonic oscillator is coupled and 7 ^ has 
the dimension of frequency. 

The steady states of the two systems are (see Appendix 

[a 

= (1-s)|0)^(0U + s|1)a(1U, 

Pst = (53) 

n + 1 \n -\-l J 

n—O ^ 


We assume that the two systems interact through the 
Hamiltonian 

Hj = na, (b^ + S), (54) 

and the evolution of the whole system is given by the 
master equation 

^ =-i[Hi,g]+C^g +C^g. (55) 

A similar type of optical master equation has been sug¬ 
gested to describe the properties of defect centers in crys¬ 
talline structures [l^. The eigenvectors of = Uz are 
| 0 )a and | 1 )^ and therefore we can construct the eigen¬ 
vectors of the superoperator C^p = [&z, p\'- 

= |0)a( 1U, d'z=-i,y=o = |1 )a(0U, 
(T;=o,i.=o = |0 )a(0U, (Tz=o,!.=i = |1)a(1U- (56) 

In the case when 1 = ±1 the Liouville operator C,"^p = 
p-\-iojA[^z^ p] has only one eigenvalue rji= ±1 = -7^/2. 
It follows from the results of Sec. uni that the steady 
state of the two-level system remains independent of the 
interaction strength H. However, for this particular case 
we are going to show the independence by a direct calcu¬ 
lation. Therefore, in order to investigate the interaction 
of the two systems and also their independent damping 
mechanism we now consider the density matrix in the 
following representation 


Q = 


Pop{a,a*) Pop{a,a*) 
Pip{a,a*) Pip{a,a*) 


\a)B{a\Bd'^a, (57) 


where (Pa = (iRe{Qf} dlmja} and Pij{a,a*){i,j € 
are quasi-distributions of Glauber-Sudarshan type 
. We are going to focus only on the off-diagonal 
element Pq,!(«;«*) and using the following identities 



b^ 


\a)i 


A 

da 




b\a)B = 




(58) 

the equation of motion (ISSl) for the representation (1571) 
yields 


dPoA 

dt 


(^iwA — dD[PoA] 


where 


D[Po,i] = -ii^B 


(A¬ 


0,1 


(60) 



da* 


da 


+ 2n 


d^ 1 
dada *. 


Po.i. 


In order to obtain the above differential equations an 
integration by parts was taken in the variables a and a*. 
The structure of Eq. (l59ll suggests a Gaussian ansatz of 
the form 


_ g-a(t)-|-6(t)a+c(t)a*-d(t)|a|^^ ^0^) 

Substituting Eq. (1611) into the Eq. (l59l) . and a compari¬ 
son of the terms proportional to I, a, a* and |ap leads 
to the following set of differential equations: 

d = —2iu)A — i^{c-\-h) — (be — d) — + 7^/2, 

i) = ioJBb — ifl (d + 2) + /2b — ^^nbd, 

c = —iwsc — in (d-I-2)-I-7'®/2 c — 7 ®hcd, 
d = 7'®d —7^hd^. (62) 

The solution of these equations is very complicated and 
is subject to the initial conditions. However, we are in¬ 
terested in the limit t —>■ oo and b(t), c(t) and d(t) tend 
to a finite limit, but the normalization of Eq. dSl) is 
characterized by a{t) « 'y^t/2. Therefore in the long 
time limit Po.i(t —>■ oo) = 0 and the decay to this limit is 
upper bounded by e”'*' The same arguments apply 
to Pi.o(t). It is immediate that the steady state of the 
spin system remains invariant 

TTB{g(t oo)} = (I - s)|0)a(0U + s|1)a(1U- (63) 


In order to gain some insight into the evolution of the 
diagonal elements, we study the average excitation of the 
harmonic oscillator (Pb) = Tr{6t5p(t)} and the popu¬ 
lation inversion of the spin system (dz) = Tr{i3'ze(i)}- 
Making use of Eq. (1551) we can derive the following equa¬ 
tions of motion for (Pb) and for (dz) 

^(6^6) = —2HIm[(dz5)] — 7^(5^5)-b 7^h 

^(dz5) = ^-iujB {dzb) - 
-b 7 ^ (2s - I) ( 6 ), 

^(b) = (S) - fll(dz), 

^(dz) = -7^(dz)+7'^(2 s - 1), (64) 

where Im[z] denotes the imaginary part of the complex 
number z. 

The obtained equations form a closed set of differential 
equations and they can be solved for a set of given initial 
conditions. However, we are interested in the behaviour 








of these quantities in the steady state Qst- Setting the 
left hand side of Eq. (1641) to zero we find the following 
solutions to the population inversion 

{^z)st = Trlo-^Pst} = 2s - 1, (65) 


and to the average excitation of the harmonic oscillator 


{b'^b)st = n + (2s - 1) 
27^ 


+4(1 — s)s- 


'7 


(7^)^+4a;| 
^ 4112 


(27^ + 'y^y + 4w| 


. ( 66 ) 


The average excitation of the harmonic oscillator (6m)st 
shows an increasing dependency on the interaction 
strength H. The most peculiar scenario occurs when 
s = h = 0, in which case the non-interacting systems 
relax to their ground state pl^ = |0)i(0|i {i G 
While this excitation extracting mechanism dominates 
both systems, due to the interaction Hamiltonian in Eq. 
(15^ the average excitation of the harmonic oscillator for 
n > 0 takes the value 


{b^b)st 


4H2 

+ 4a;|' 


(67) 



n// 


FIG. 2: The average excitation of the harmonic oscillator in 
the steady state as a function of We set 7 ^/ 7 ^ = 1, 

h = 0 and s = 1/2. The solid line with = 1 and the 

dashed line with 01 ^/ 7 '^ = 5 shows an increase in the average 
excitation, but this increase is mediated by the frequency of 
the harmonic oscillator . 


Fig. [2] shows that the average excitation of the har¬ 
monic oscillator starts at the value n and increases as 
a function of H. This increase depends on the relax¬ 
ation parameters of the spin system and of the har¬ 
monic oscillator. The found behaviour suggests that 
the model can only describe physical scenarios whenever 

We also investigated numerically the decaying mecha¬ 
nism of the whole system. We considered the initial con¬ 
dition to be |0 )b- Thereby choosing h = 0 and 



0 1 2 3 4 5 6 

/t 

FIG. 3: The norm distance between the evolved state p^(t) = 
TrB{^(t)} and the steady state of system T as a function of 
7 ^t. We set = 1, h = 0, s = 1/2 and = 

= 10. The solid line with H/y'^ = 5 and the dashed 
line with fl/y'^ = 0 shows the difference between coupled and 
uncoupled system. 


s = 1/2 the only elements affected by the dynamics are 
the off-diagonal elements A(0|e(t)|l)A and A(l|p(t)|0)A- 
We simulated the harmonic oscillator with ten basis vec¬ 
tors which are enough because the initial condition is 
the steady state for system B and the chosen coupling 
strength O is not big enough to drive the harmonic os¬ 
cillator towards higher excitations than 10 and the con¬ 
tributions from these states are negligible. Fig. [3] shows 
that the decay e“^ of the uncoupled system A, i.e., 
H = 0, is the upper bound for the off-diagonal elements 
which is in agreement with our general findings in Sec. 

EH 


C. The optomechanical system 


In this section we study more in detail the model which 
motivated our work, see Sec. m The system is composed 
of two harmonic oscillators, i.e., two infinite systems, and 
a direct approach to find out the steady state of the joint 
system requires the solutions of an infinite number of 
recurrence relations. Therefore our result in Sec. EH is 
useful to identify some properties of the steady state. The 
Hilbert spaces are Ha = rs(C) and Hb = rs(C) with 
orthonormal basis {|n)A}nGN for system A (optical mode) 
and {|n.)s}nGN for system B (mechanical mode). The 
uncoupled systems are governed by the master equation 


dt ’ dt 




( 68 ) 


with the Liouville operators defined in Eq. @. We re¬ 
call that a (a^) and b (6^)are the annihilation (creation) 
operators of the quantum harmonic oscillators A and B 
with frequency uj and frequency v respectively. The coef¬ 
ficients with dimensions of frequency characterizing the 
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dissipative dynamics are k and 7 for the optical and 
mechanical mode respectively. The temperature of the 
photonic reservoir is proportional to n and to m in the 
phononic case. The steady states of the separate systems 
are (see Appendix |B| 

''A ^ 

n—0 

-| 00 y _ X n 

A® = ;^E(5rrT) 

n—0 ^ 

The two systems interact with each other through the 
optomechanical interaction Hj = ga^a and the 

evolution of the whole system is given by the master 
equation ©• 

Connecting with the formalism in Sec. IIIII we note 
that in this case a and for all n £ N |n)yi is an 

eigenvector with eigenvalue n. The eigenvectors of the 
superoperator C^p = [a'^'a,/)] are 


\n)A{n\ 


— 


\n + l)A{n\A, I > 0 

\n)A{n + 1 \a, I < 0 


n G N. 


(70) 


In the subspaces defined by the eigenvectors I 7^ 0 the 
smallest eigenvalue is pi = —k\1\/2 (see Aponedix lAl) . 
Applying the results of Sec. |nT]we get 

Trsigst} = ptf (71) 


In this particular example it is hard to work out an al¬ 
ternative proof for the above relation as we did in sec¬ 
tion |TVBj In general Tr^jpst} 7^ Pst and using the same 
approach as in the previous section we are going to an¬ 
alyze the average excitations (a^a) = Tr{aldp(t)} and 
(b^b) = Tr{b^bg{t)}. The time evolution of the density 
matrix g(t) is governed by Eq. ([2]) and the following 
equations of motion can be derived for (a^a) and (b^b) 


dt 


(St6) 


dt 


(d^db) 



dt 


((aU) 



—2glm[(a^a6)] — plpb) + 7m, 

-/£\ 7 + 2(K-bw) + 

Kn{b) -^- [a'ab) — ig{[a'a) ) 

-ii/(6) - ig{d'd) - -(6), 

Kh —2tt;((d^a) )-|-K(4h-I-l)(a^a), 

—K{d^d) -I- Kh. (72) 


The steady state solutions of Eq. (EH) are 


(a^d)st 


= n, 


= m + 


An^g^ 
7^ -I- 


(73) 

4h(h -I- 1)(2 k -I- 7 ) 5 ^ 

7 [(2k -I- 7)2 + 4 ^/ 2 ] 


The above relation shows that the average excitation of 
system B is increasing as a function of g. The only case 


when it remains the same is for n = 0, i.e., system A 
is relaxed to the ground state. So, whenever system A 
is relaxed to a finite temperature Gibbs state the aver¬ 
age excitation number of system B is increased above 
the average excitation given by its own damping mech¬ 
anism. These results resemble the same effect found in 
Sec. IIV Bl The parameter region where this model works 
is constrained by the relation v ^ gnt/ K/7. This is usu¬ 
ally fulfilled, because in experimental scenarios h « 0 



FIG. 4: The norm distance between the evolved state p^(t) = 
TrB{^(t)} and the steady state of system A as a function of 
nt. The dashed curve corresponds to the free evolution of A, 
i.e. g = 0. The full line with g = 0.9k shows the difference 
between the coupled and the uncoupled system. The initial 
state for both systems A and B is a coherent state |q}, with 
a = 0.15. The parameters are: cj = 10k, v = 1.5k, 7 = 0.9k, 
h = 0.015 and m = 0.1. 


Eig. |4]shows the relaxation of the optical mode to the 
steady state. We numerically evaluate the trace norm 
||p^(t)—p^lli to measure the distance of the evolved state 
to the steady state in Eq. inH The dashed line shows the 
free scenario where the optical mode is solely damped 
from an initial coherent state |a) with a = 0.15. The 
full line shows the reduced density matrix of the optical 
field after the optomechanical interaction. Both modes 
are taken initially in a coherent state with a = 0.15. One 
can note that the trace norm vanishes for both situations, 
meaning that the optical system always relaxes to the 
steady state . Even though the optical steady state is 
the same, the relaxation mechanism changes due to the 
interaction with the mechanical mode. 


V. CONCLUDING REMARKS 

We have discussed the steady state of two interacting 
quantum systems each of them having an independent 
non-unitary Markovian evolution. The interaction be¬ 
tween these two systems was considered in a way that 
it commutes with one of the system Hamiltonians and 
that the unitary operator generated by it leaves the Li- 
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ouvillian of the non-unitary dynamics invariant. It is 
demonstrated that the steady state of the whole system 
has an interesting property, the trace over the other sys¬ 
tem’s degrees of freedom results in a steady state for the 
commuting system which is also the steady state of this 
system without the interaction. In plain words, the com¬ 
muting system relaxes always to the same steady state, 
independently of the interaction strength. 

The method applied here is based on the knowledge of 
the spectrum of the commuting system’s Liouville oper¬ 
ator. The eigenvalue equations of these generators are 
known only in a few cases. Therefore we restricted our 
statement to the case where the system with the com¬ 
muting properties is composed of spins and harmonic os¬ 
cillators with damping mechanisms defined in Eqs. 
and though the presented method is general enough 
to deal with other systems as well. In the case when this 
system is finite (more than two levels) the solutions to 
the eigenvalue equations can be done either with analyt¬ 
ical formulas or with numerical simulations. However, 
we chose the model of a two-level (spin) system with 
decay processes due to its frequent application in var¬ 
ious physical scenarios. In the case of infinite systems 
our knowledge is restricted to the harmonic oscillator. 
The decoherence mechanism of a harmonic oscillator can 
be described by an optical master equation whose eigen¬ 
value equations are known Q. This is the reason why 
this model is our second choice. These damping mod¬ 
els with their Liouvillian generators are considered in 
the decomposition of the commuting system, although 
we allow internal interactions between the systems. For 
example: in the case of a spin and a harmonic oscilla¬ 
tor we have the two independent damping mechanisms 
and a Jaynes-Cummings type interaction between them; 
for two spins and a harmonic oscillator we have three in¬ 
dependent damping mechanisms and a Tavis-Cummings 
type interaction between them; for a spin and two har¬ 
monic oscillators we have again three independent damp¬ 
ing mechanisms and a dipole coupling in the Lamb-Dicke 
regime, an interaction used to describe trapped ions with 
center-of-mass motion inside an optical cavity. In our 
analysis we do not have to specify the other system. We 
comment here that in the case of a Caldeira-Leggett mas¬ 
ter equation [2^ the solution to the eigenvalue equation 
is still unknown. The proof uses theorems from the the¬ 
ory of one-parameter dynamical semigroups. The reason 
to apply these results is the Liouville operator of the har¬ 
monic oscillator which is an unbounded operator acting 
on the Banach spaces of trace-class operators. We use 
the fact that dynamical semigroups are contractive semi¬ 
groups on the subspace of trace-class operators whose 
elements are self-adjoint operators. We slice up this sub¬ 
space into further closed subspaces and derive in each 
of them an inequality for the norm of the evolved den¬ 
sity matrix with the help of the Trotter product formula. 
These inequalities allow us to identify which elements of 
the density matrix tend to zero in the limit of inhnite 
time. The projections which splice up the subspace are 


defined by those vectors which are eigenvectors of the 
interaction’s superoperator. 

We considered three examples of this scenario: a spin- 
spin, a spin-harmonic oscillator, and a two harmonic os¬ 
cillator system. In the first example every detail about 
the time evolution and the steady state can be worked 
out analytically due to the small dimension of the Hilbert 
space. For the other two systems the situation is more 
complicated and in order to understand the evolution of 
the system we derived closed differential equations for 
the averaged excitations (harmonic oscillator) and the 
population inversion (spin). We focused on the steady 
state solutions which show that the interaction increases 
the average excitation number only for the system with¬ 
out the commuting properties. In the case of the op¬ 
tomechanical model, it turns out that this effect is unde¬ 
tectable in current experimental realizations due to the 
weak coupling between the optomechanical resonator and 
the single mode of the radiation field and the low value 
of thermal photons in optical frequencies. 

Finally let us make some comments on the limitations 
and possible extensions of our approach. The results of 
this work are based on a class of coupled systems with 
independent Markovian master equations, models which 
can explain certain physical scenarios like the actual op¬ 
tomechanical experiments. However, it is known that 
strong interaction between two coupled dissipating sub¬ 
systems leads to serious defects, whenever independently 
derived non-unitary evolutions are added to the unitary 
evolution [27l - l35j |. If these defects are present in the so¬ 
lutions, then one must start the theoretical aproach on 
a microscopic level and derive a new dissipative model. 
The structure of these master equations depends on the 
interaction between the two subsystems. For example, in 
Ref. [13 the two subsystems are harmonic oscillators and 
their creation and annihilation operators are mixed up in 
the generator of the master equation, which leads to an 
open problem of determining the spectrum of this Liou¬ 
villian. Our statements are based on the requirement of 
independent master equations in order to keep the non¬ 
commuting system arbitrary, however when the details of 
this system are known and they fulhll the assumptions of 
Sec. nni then the steady state of the commuting system 
remains invariant, see for example the strongly coupled 
dissipative optomechanical model [l^ . 

As a final word we think that the presented analysis 
may offer interesting perspectives for actual and future 
models with Markovian master equations. 
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Appendix A: Eigensystem of the Liouville operator 

In this appendix we present the eigensystem of a 
damped two-level system and a damped harmonic os¬ 
cillator. Due to the fact that the Hilbert space of a har¬ 
monic oscillator is infinite dimensional we are going to 
summarize some basic facts about the adjoint of a Ba¬ 
nach space operator, a notion needed for the definition 
of the eigenvectors. 


1. Eigensystem of Us 

In the main text we considered the following master 
equation for a spin system 

Csp = - m) {&-pa+- \{a+a-,p}) 

-bym ((T+ - i {ct_ CT+, /3}) . (A1) 

The corresponding Liouville operator is non-Hermitian, 
but it can still be diagonalized leading to the dual eigen¬ 
value equation 


^SPi — ^iPi: ^sPi — ^iPi- 

The elements pi are the four right eigenvectors given by 

po = \i+ Pz = a-z, p± = a±, (A3) 

with the corresponding eigenvalues 

Ao = 0, , A. = -7, A± = 7/2. (A4) 

The adjoint of the Liouville operator is 

= 7(1 - m) {o+P^- - i P}) 

-b7TO {a+pa- - i {a-a+,p}) , (A5) 

and its eigenvectors are the following four elements 

P0=i, Pz = P±=CT±. (A6) 

The eigenvectors of £s and £g belong to the space of 
2x2 matrices Al 2 (C) and form a biorthogonal system 
as it can be seen that they fulfill the relation 

(PiJj) = Tr|/5jPi} = kj (A7) 

with the Kronecker delta for all i,j G {0, z, -b, —}. 


2. Eigensystem of £ho 

The master equation describing a Markovian damped 
harmonic oscillator can be written in the following form 

-Chop = lin + I) (apa'^ - p}) 

+jn pd — ^{aa^, p}) . (AS) 


The eigenvalue equations of £ho and £{jq read 

^HOPrijk — ^n,kPn,k-; ^\iQPn,k — ^n,kPn,k- (A9) 


The first equation is solved by the following right eigen¬ 
vectors (see Ref. i) 


. I l>c|+fc 

Pn,k = a' 2 


(-1)" . r(|fc|) p- 

(n-|-l)'=+i ■ \fi+l J ^ 


n+1 


|fc|—fc 
2 


(AlO) 

where is a generalized Laguerre polynomial and 
: /(did) : is the normal ordering of /, a function of the 
number operator d^d. The eigenvalues are 


An,fc — 7 



n G No, fc G Z. (All) 


£|jq is a representation of the adjoint of £ho with the 
help of the trace as a linear functional and can be found 
to be 


= l{n + 1 ) (dlpd - i{dld, p}) 

{dpd} — i{ddl,/5}) , (A12) 

and its eigenvectors are given by 

(A13) 


It can be verified that they are orthogonal with respect 
to the inner product 


{Pn^ki Pn' ,k'} — j — ^n,n' ^k,k': 

and this is consistent with 

{^}loPn,k-j Pn,k) — {Pn,k: ^liOPn,k) — ^n,k- 


(A14) 


(A15) 


3. Comments on the eigensystem of £ho 

Banach space structure.- In the case of the damped 
harmonic oscillator one has to deal with an inhnite di¬ 
mensional Hilbert space. The Liouville operator is a map 
on the set of trace class operators which do not form a 
Hilbert space, but do form a Banach space. Geomet¬ 
rical problems are inherent in the concept of a Banach 
space and therefore a slightly different notion of adjoint 
operators is needed. 

Let E, F be Banach spaces and if T : i? —>■ F is a 
bounded linear operator then we define its adjoint by 

T' :F' ^ E', T'{ip) ^ipoT,ip€F' (A16) 

where T' is mapping bounded linear functionals on F to 
bounded linear functionals on E and E'{F') is the dual 
of E{F). The set of trace class operators on a Hilbert 
space H is defined as 

Bi{n) := {x G B{n) : Tr 


x^x [• < oo}. (A17) 








12 


If £ : —)■ is a bounded linear operator then 

the adjoint C is defined on the dual Biiji)'. Due to the 
fact that the set of trace class operator form a two-sided 
ideal in (see Theorem VI. 19 in Ref. [1^) we have 

an isometric isomorphism 

Bi{ny'^B{n). (ais) 


and therefore they are not in "H. In the case of X 2 , it 
can be noted that ax 2 £ V. but a^aa::2 ^ H. With this 
observation we can conclude that Dom(a’'’a) C Dom(a) C 

n. 

In the case of the Liouville operator and for the sake of 
simplicity, let us consider the case of zero a temperature 
bath (n = 0). In this case the Liouvillian reads 


Thus, we may translate the adjoint C' to an operator 
: B{'H) B{'H) satisfying 

TT{{C^)^p} = Ti-{yHCp)} (A19) 

for all y G S(H) and p G Bi{'H). 

We can already conclude that the eigenvectors of 
can be understood as the linear functional (p for a fixed 

Pn,k 


^[Pn,k]{p) = Vp G Dom(/:Ho)- (A 20 ) 

Domain.- When dealing with infinite dimensional sys¬ 
tems, the domain of the operators is a non-trivial con¬ 
cept that has to be handled carefully as we are going 
to elucidate in the following lines. The Hilbert space 
of the harmonic oscillator is the symmetric Fock space 
Ti, = rs(C). Every element can be expressed in terms of 
the orthonormal basis {|n)}„gN as 

OO OO 

a: = ^ a;„|n) G "H, ^ < oo. (A 21 ) 

n—0 n—0 

The creation and the annihilation operators (a^, a) have 
the domain 


Dom(a^) = Dom(a) = |a; : n|a:„p < oo| (A22) 

n—1 

which is dense in T-L. The number operator a^a has a 
different and more restricted domain given by 

OO 

Dom(a^a) = jx : < oo|. (A23) 

n—1 

In order to exemplify this difference, let us define the 
following three vectors in "H: 



= (A24) 

n=0 

These vectors have the following properties 

Xi ^ Dom(a), xi ^ Dom(a^a), 

X2 G Dom(a), X2 ^ Dom(a^a), 

X3 G Dom(a), x^ G Dom(a^a). (A25) 

This is because even when x\ is a normalized vector, 
axi ^ Ti, and a^axi ^ Ti, are not normalized anymore 


Cop = apa) — - {ayap + pa^a) . (A26) 


In order to show that not all elements of Bi{'H) are in the 
domain of Cq , we consider the following density matrix 


. 6 ^ 1 I \/ I 

n,m—l 

(A27) 

The normalization condition is fulfilled as = 

tt^/6. Applying Cq one gets 


P* 




TT" ^ V«+T 

n—0 


(|n)(0| + |0)(n|). 


E 


^ n7^=i \ ^/{n + l){'m + l) 


n -\- m 
2nm 


|n)(m|. 


Now, let us remember that the normed space of Hilbert- 
Schmidt operators is a subset of all bounded linear oper¬ 
ators defined on a Hilbert space 'H and at the same time 
it contains the set trace class operators, i.e. 


Bi{n) c B2{n) c B{n). (A28) 


In order to determine if p* is a Hilbert-Schmidt operator 
we evaluate 


Tr{ptp*} = E 


;^i V\/(’^ + l)("i+ 1) 


n m 
2nm 


+ 2E- + I 


(A29) 


which is clearly diverging. Therefore p* ^ B 2 {'H) im¬ 
plying that p* ^ Bi{'H). So £0 in Eq- (IA26D cannot be 
defined on the whole space Bi{'H), but only on its domain 
Dom(£o). The properties of Dom(£Ho) are not studied 
here, but one can at least state that 

Dom(£Ho) CHi(H), 

which implies that the bounded linear functionals on 
Dom(£Ho) regarding to the norm of Bi{'H) are more nu¬ 
merous than the bounded linear functionals on Si('H). 

Linear independence and completeness.- The eigen¬ 
vectors pn,k due to the normal ordered form : 
Ln\zd^d)e~^°'^°' : with z G (0,1] are trace class oper¬ 
ators. Nonetheless, {a)a)™'pn,k is also a trace class op¬ 
erator for any finite value of to. All the eigenvalues of 
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£ho are distinct which implies that the eigenvectors pn,k 
are linearly independent. We have a countable infinite 
set of eigenvectors whose linear independence property 
is equivalent to show that all possible finite selection of 
eigenvectors are linearly independent. In order to show 
this, let {/3i(i),/ 3 j( 2 ),..., Pi(fe)} be a finite set of eigenvec¬ 
tors with corresponding eigenvalues Ai( 2)7 Ai(fc)} 

and let ci,C 2 , ...,Cfc be scalars such that 

ClPi(l) + C2Pi{2) + + CkPi{k) = 0. (A30) 


of a frequency. To ensure that Eq. m describes the 
evolution of a density matrix we require the condition 
72 < 71- 

To investigate in more detail the steady state of the 
damping, we consider the density matrix of the field in 
number state representation \n) (n £ No) 

OO 

Pst = ^ Pn,m\n){m\. (B2) 

n,m=0 


£ho has non-degenerate eigenvalues and applying 
■A^(£ho) = (£ho - Aq2))(£HO - Ai(3))...(£HO - K(k)) 
to the above equation results that all terms except the 
first are annihilated, resulting 

ciiV(A,(i))ft(i) = 0, iV(A^^i)) ^ 0, (A31) 

and since pi(i) ^ 0 it follows that ci = 0. Similarly, 
C 2 = ... = Cfc = 0 , which proves the linear independence 
of {Pi(i)j/ 5 i( 2 )j ■•■7 Pi(fc)}- This proof can be applied to all 
possible finite sets of eigenvectors. The eigenvectors are 
not orthogonal in regard to the Hilbert-Schmidt inner 
product. The set of eigenvectors has maximal possible 
cardinality, so {pn,k}n&i,k& is a basis. 

Spectrum.- Finite rank operators are dense in the set 
of trace class operators (see the Corollary after Theorem 
VI.21 in Ref. (Sg). The eigenvectors pn,k are finite rank 
operators and they form a basis which implies that £ho 
is a closed densely defined operator. The point spec¬ 
trum of £|jq coincides with the point spectrum of £ho 
and this means that £ho has no residual spectrum (see 
Proposition 1.12 in Chap. IV of Ref. [111). Now we 
know the point spectrum and the residual spectrum, but 
we have not yet discussed all possible approximate eigen¬ 
values. A is an approximate eigenvalue if there exists a 
sequence {in}nGN C Dom(£Ho) such that ||in||i = 1 
and lim„_>oo ||£Hoin — Ai„||i = 0. The latter equation 
can be expanded in regard to the basis {pn.fe}nGN,feGZ and 
yields 

lim II (A„,fe - A) Tv{gn,kin}pn,kh = O 7 

n—¥oo ^ 

n,k 

which is fulfilled only when A„,fc = A, meaning that £ho 
has only a point spectrum. 


Appendix B: Steady state of a damped harmonic 
oscillator 


We set the left-hand side of Eq. m to zero, and taking 
matrix elements of the density matrix we get 

^1 ! - 

0 = - -^(m-|-n)p„_^-|-7iv'(n-|-l)(m-|-l)p„+i,m+i 

'Y 2 

- (m-I- 1 -I-n-I- + 72t/«TOPn-l.m-l. 

(B3) 

The matrix elements are coupled via a second order re¬ 
currence relation to their nearest neighbours in each of 
the diagonals of the density matrix. 

First, we consider the main diagonal (m = n) and Eq. 
(IB3I) simplifies to 

_ 7in-I-72(n-I-1) 72n 

- yi(n + l) + 

(B4) 

The solution can be determined by induction, thus yield¬ 
ing 


Pn.n = ( 72 / 71 )" Po ,07 11 > 0, (B5) 

The normalization condition Trjp} = 1 results in 



With the geometric series 


00 / \ n 

=_ 1 _, 

1 - 72 / 71 ’ 


(B6) 


(B7) 


one finds that po_o = (71 — I 2 )Hi- 

Now, we consider only the upper diagonals m > n 
due to the self-adjoint property of the density matrix ps- 
Defining I = m — n> 1 and 72 = £71 with e £ [0,1) Eq. 
(|B3p takes the form 


Let us consider the model of a damped harmonic os¬ 
cillator given by the master equation 


dp 

dt 


71 

2 

72 
2 


(a)dp — 2a/5al -I- pa^d') 
(afilp — 2al/5a -|- pda)) 


(Bl) 


where d (&!) is the annihilation (creation) operator of 
the harmonic oscillator. 71 and 72 have the dimension 


Pn,n+l —l’;,nP0,Z7 

Hn+l = (/" + Hn - <^Krln-i (B8) 
where we introduced the following abbreviations 
/" = (n + 1/2) /H[n l)(n -|- Z -|- 1), 

= (n -b 1 -f ^2) /\/(ii + l)(ii + ^ + 1)7 
h) — H 11(11 T 0 /\/(ii + l)(ii + ^ + 1 ), (B 9 ) 
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which are all positive and fulfill the following inequalities 

fP-hf>0, g->l, Vn>0,/>!. (BIO) 

For the first few terms it can be shown that they have 
solutions of the form 


Thus, we have found that all the solutions to Eq. (IB3I) 
for / ^ 0 fulfill the inequality 


I y ] Pn,n+l I > \P0,l I 

n=0 



V 7re^/^(n + 1) 


(B19) 


n 

(Bll) 

i=0 

with coefficients obeying the properties 

a"’* > 0, Vn > f > 0. (B12) 

We apply a proof by induction to show that takes 

the form of Eq. (IBlip with the conditions of Eq. (IB12P . 
First we use Eq. (|B 8 p which results in 


<n+i = E +E5r«r’'e*+' - E 

2=0 2=0 2 = 0 

(B13) 

and translates to equations for the coefficients: 

= fraf, (B14) 


Now, applying the conditions of Eqs. (IBlOp . and (|B12ll 
we find that 


n+l,Z 


n,Z 


n+l,Z 

^i+l 


(B15) 


and it follows immediately that every ^ can be cast in 
the form of Eq. (IBllI) . 

As all coefficients are positive, the sum is always larger 
than the zeroth order coefficient 


n—1 


r 


0 _ 

l,n 


If/; 


r('i + l) / i! 

r(|) y n\{n + l)V 


(B16) 


which we have expressed in terms of the gamma function 
r( 2 ;), and multiplied by po.i solves the first order recur¬ 
rence relation obtained from the Eq. (IB3I) in the case 
when e = 0 which implies 72 = 0 . 

These considerations allow us to set a lower bound to 
the following sum: 

00 00 00 

\'^Pn,n+l\ > \P0,l\J2^ln > IPO.zI (B17) 

n—0 n—0 n—0 

as the relation > r)* „ can be verified by noting that 
each term in the product of Eq. (IB16I) has the property 
/r > fi which can be checked from Eq. (IB9I) . Stirling’s 
formula in an inequality form [ 3 ^ [38l | 


The right-hand side of the above inequality is divergent 
for |po,z| > 0 because it contains a harmonic series. 

In the next step we show that Eq. (IB19I) is upper 
bounded by 1 due to some general facts of operator the¬ 
ory. Let us consider an arbitrary separable Hilbert space 
TL with an orthonormal basis (e^). The set of Hilbert- 
Schmidt operators B 2 {'H) with the Hilbert-Schmidt inner 
product 

{X,Y) = E(^e„fei) = Tr{lty}^ VX,f G 02(22) 

(B20) 

is a Hilbert space and it is also a two-sided ideal in 
B{'H){see Theorem VI.22 in Ref. [3^1: 

PVII 2 , ||Vi ||2 < 00 , V G 02(22), A G 0(22). (B21) 

In our case the Hilbert space 22 = rs(C) and any den¬ 
sity matrix p G 02 ( 22 ) (||/ 3||2 < 1) and G 02 ( 22 ) 
( 11 / 5^/^112 = 1 ). 

Now let us consider the right shift operator 

p|n) = |n -I- I), n > 0, 

Sr G 0(22). (B22) 

Any Ah power of the right shift operator 0)^ is also a 
bounded operator. Then 

00 

|Tr(S!,ft)| = I Pn,n-\-l I (B23) 

n=0 

and 0(jPs G 02 ( 22 ). The Cauchy-Schwarz-Bunyakovski 
inequality for the Hilbert-Schmidt inner product gives: 

|Tr(XtF)| = |(V,F)| < ||V||2||y||2, VX,y G 02(22). 

(B24) 

We observe that for the Hilbert-Schmidt operators X = 
pIP and Y = 

|Tr(4p,)| < \\pII%\\S^rpII% = 1, (B25) 

where we used the relation 0 )^ = 2 , the identity 

operator on 22. Therefore, Eq. (IB19I) obviously leads to 


00 1^1 

1 > IE ^"."+'1 > lpo,/|+^=|po,/l E ~rT- 

„=o V7rei/3 


_ / jl\'n _ / Tl\'^ 1 

■v2TTn y—j < nl < v27rn y—j (B18) 

and the relation r(n -|- 1/2) = Y^(2n)!/(4"n!) yield 

g e 6" 1 

y/Trn{n + 1 ) V'7Te^/^{n + 1 ) 


which is true if and only if po^i = 0. If po^i = 0, then 
all the upper diagonal elements of pst are zero. It is 
immediate from the self-adjoint property of pst that all 
lower diagonal elements are also equal to zero. Finally, 
we can conclude that pst has only non-zero elements in 
the diagonal which are given in Eq. (IB5I) . 
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